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We study the quenching dynamics of a many-body system in one dimension described by a 
Hamiltonian that has spatial periodicity. Specifically, we consider a spin- 1/2 chain with equal 
XX and yy couplings and subject to a periodically varying magnetic field in the z direction or, 
equivalently, a tight-binding model of spinless fermions with a periodic local chemical potential, 
having period 2g, where q is a natural number. For a linear quench of the magnetic field strength 
(or potential strength) at rate 1/r across a quantum critical point, we find that the density of 
defects thereby produced scales as 1/r'''^''''^-', deviating from the \l\pr scaling that is ubiquitous 
to a range of systems. We analyze this behavior by mapping the low-energy physics of the system 
to a set of fermionic two-level systems labeled by the lattice momentum k undergoing a non-linear 
quench as well as by performing numerical simulations. We also find that if the magnetic field is 
a superposition of different periods, the power law depends only on the smallest period for very 
large values of r although it may exhibit a cross-over at intermediate values of r. Finally, for the 
case where a zz coupling is also present in the spin chain, or equivalently, where interactions are 
present in the fermionic system, we argue that the power associated with the scaling law depends 
on a combination of q and interaction strength. 

PACS numbers: 64.70. Tg, 75.10.Jm, 71.10.Pm 



I. INTRODUCTION 



The quenching dynamics of quantum many-body sys- 
tems has become a topic of active research in recent 
yearsir— . In particular, there has been a focus on the 
effects of slow quenching across a quantum critical point 
(QCP). Here a parameter A in the Hamiltonian describ- 
ing a quantum system is varied linearly (for instance) 
in time at a rate given by l/r, such that the system 
starts in the ground state far away on one side of a QCP, 
crosses the QCP at a critical parameter value Ac, and 
ends far away on the other side of the QCP. The final 
state of the system is in fact dominated by the regime 
close to the QCP—"— where the correlation length of the 
system diverges as |A — Ad"'^ and the relaxation time (or 
the inverse of the energy gap between the ground state 
and the first excited state) diverges as |A — \(\~^^ . No 
matter how small the quenching rate 1/t, the relaxation 
time for the system is larger than r because at the QCP, 
there exist modes with energies arbitrarily close to zero. 
The quenching process near Ac is not adiabatic for such 
modes, and the system does not reach the ground state 
of the final Hamiltonian. The final state contains a finite 
density of defects which scales as a power of r for small 
l/r given by the Kibblc-Zurek (KZ) form 
The power thus depends on only three quantities; the 
spatial dimensionality d, and critical exponents v and 
z. For a large class of translationally invariant systems, 
the Hamiltonian decouples into several pairs of momen- 
tum modes of opposite sign. In such systems, quench- 
ing dynamics and power-law scaling can be analyzed by 
mapping the problem to the dynamics of the Landau- 
Zener transition of a two-level system dictated by a time- 



dependent tuning parameter—. In several such instances 
in one dimension, one has z = v = \, giving rise to the 
ubiquitous 1/t^/^ scaling of the density of defects scales. 

While the above scaling law holds for a range of sys- 
tems, some generalizations and deviations thereof are 
coming to light. For instance, if the Hamiltonian is var- 
ied across the QCP non-linearly in time, as (t/r)", the 
scaling law becomes i/r'''^"/(^'^"+i)^ao. Effectively, the 
non-linear quench modifies the exponent v to vol. As an- 
other instance, it has also been shown that the power 
law can depend on certain topological features of the 
systeroi^iiiiil. In particular, in Refs. [13 and [13, some of 
the authors of the present work have studied quenching 
dynamics in a two-legged ladder version of the Kitaev 
model. As with the parent Kitaev model on the hexag- 
onal latticei^"— , the two-legged ladder has a large num- 
ber of sectors (growing exponentially with the number 
of sites) which are distinguished by the eigenvalues (±1) 
of a number of local Z2-valued conserved quantitie o^^'^^ . 
The ladder was mapped to a fermionic p-wave supercon- 
ducting system which has recently emerged as an exciting 
topological system in its own righti^. It was shown that 
while quenching through QCPs in certain sectors yields 
the standard l/r^/^ scaling of the defect density (also re- 
lated to the total probability of excitations) , certain pe- 
riodic patterns in the Z2 invariants give rise to a I/t'^^^ 
scaling. It was further conjectured in Ref. [13 that more 
general scaling laws of the form \lT''/ii+^) for the slew of 
integers q = 3,4, •• • may appear in some other sectors, 
but explicit examples of such sectors were not found in 
that model. 

In this work, we show that the presence of a period- 
ically varying parameter in an otherwise spatially ho- 



mogencous system provides an excellent route for new 
and interesting power-law quenching behavior, inclusive 
of the l/r''/(^+^) scaling. We demonstrate the role of 
periodicity in the specific case of a one-dimensional spin- 
1/2 lattice model having homogeneous and equal nearest- 
neighbor XX and yy couplings and subject to a trans- 
verse periodically varying magnetic field having the value 
hcos{Qn + (p) at site n. By the Jordan- Wigner trans- 
formation, we equivalently study a tight-binding system 
of non-interacting spinless fermions in one dimension in 
which the local chemical potential varies as h cos{Qn+(t)). 
For Q = Tv/q (where q = 1,2, ■ ■ ■), the system exhibits a 
spatial periodicity of 2q, fragmenting the Brillouin zone 
into 2q regions that are coupled to one another by the 
periodic potential. In contrast to the aforementioned 
case of pairwise coupling of modes of opposite momen- 
tum, the periodicity presents a more complex structure in 
which a quench can cause probability amplitudes to shift 
between these 2q modes. While the total post-quench 
excitation probability is still dominated by the QCP at 
h = and the gapless points at momenta k = ±7r/2, it 
depends on the intermediate paths taken by the matrix 
elements connecting the various fragments of the Bril- 
louin zone; the parameter controls the relative phase 
between the different paths. Using perturbation the- 
ory to q-th order for small h, we explicitly illustrate 
this point for a linear quench where the amplitude h 
varies linearly in time as t/r so as to go from — oo to 
oo. The perturbation provides an effective low-energy 
Hamiltonian for momentum modes around k = ±7r/2 
and reduces the dynamics to the generalized Landau- 
Zener evolution of a set of two-level systems in the pres- 
ence of a time variation, reminiscent of the non- 
linear quench. The behavior of pk, the excitation 
probability for each momentum mode, becomes increas- 
ingly complex as the value of q increases; however, a 
simple scaling analysis shows that the total excitation 
probability only gets a contributions from modes lying 
within a range of k given by Ak ^ (cos((7(/))/r'')^/'''+^) 
or {sm{q(l))/T''y^^'^~^^\ depending on whether q is odd or 
even, respectively. Hence the total excitation probability 
for the quench yields the desired l/r^/^*+^^ scaling law, 
multiplied by (cos(#))i/(«+i) or (sin(g(/)))i/(«+i) . We 
corroborate and expand our analytical arguments using 
numerics. It turns out that one can think of the effec- 
tive low-energy theory as either describing a system with 
V = 1 and a Hamiltonian varying in time with a power 
a = g, or as a system with v = q with a Hamiltonian 
varying linearly in time (a = 1). In all cases, we obtain 
a KZ power law of the form where d=l, 
z = 1, and va ~ q.Our results thus show that internal 
mode structure, when combined with critical behavior 
and quenching, can give rise to rich out-of-equilibrium 
dynamics and scaling. 

In another line of investigation, it has been shown 
that the simple post-quench 1 / y'r scaling can have dras- 
tic modifications due to a completely different reason, 
namely, interaction a^^'^^ . In Ref. [H, an interacting sys- 



tem of bosons in the presence of a periodic potential per- 
mitting one boson per potential minimum was analyzed 
within a sine-Gordon framework. It was shown that a 
quench in the strength of the periodic potential results 
in the density of defects having a power law scaling with 
an exponent that depends explicitly on the interaction 
strength. In the context of our work, interactions in a 
sense are an extreme limit of coupling between momen- 
tum modes. Borrowing from the analysis of Ref. [l2l . 
we extend our studies of periodic structures to include 
interactions in our fermionic system, or equivalently, to 
include a zz coupling in the spin chain system. We argue 
that the scaling exponent for the total excitation proba- 
bility now depends on both interactions and the periodic- 
ity q and that this result is also valid for a generalization 
of the studies in Ref. [H to the case of q bosons per 
potential minimum. 

We present our studies as follows. In Sec. II, we in- 
troduce the spin and fermion models, express the peri- 
odic term in the basis of Brillouin zone modes and derive 
the effective low energy Hamiltonian using perturbation 
theory. In Sec. HI, we analytically derive the scaling be- 
havior of the system for a linear quench in h and present 
numerical results of quenching simulations for different 
values of q, phase ^ and quenching rate l/r. To un- 
derstand the detailed dependence of on k, we study 
in the Appendix a two-level problem in which the time- 
dependent term in the Hamiltonian as \t\'^sgn(i) (where 
sgn denotes the signum function). In Sec. IV, we study 
more complex behavior of the periodic potential such as 
a superposition of two commensurate periodic functions; 
we find that the power law scaling for the total excita- 
tion probability is governed by the smallest period for 
very large values of t although the power law may show 
a cross-over at intermediate values of r. In Sec. V we 
consider the effect of interactions and we conclude with 
general remarks in Sec. VI. 

II. THE MODEL AND ITS ENERGY 
SPECTRUM 

Our starting point is a one-dimensional spin-1/2 model 
placed in a transverse magnetic field whose Hamiltonian 
is given by 

N 

H = - J + «+i + (1) 

n=-N 

where cr^ [a — x, y, z) denote the Pauli matrices at site n, 
and we are eventually interested in the thermodynamic 
limit — > oo. (We generally set h, the exchange coupling 
J and the lattice spacing a equal to unity. When we 
introduce a quenching time r, large or small values of t 
are as compared to h/ J.) Note that commutes 
with the Hamiltonian. This system can be mapped to 
a model of spinless fermions using the Jordan- Wigner 
transformation^^. At any site n, we map a spin state 
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with crfj = 1 or cr^ = — 1 to the presence or absence of a 
spinless fermion at that site; this is done by introducing a 
fermion annihilation operator c„ at each site, and writing 
the spin at that site as 

< = 2cicn - 1 = 2p„ - 1 

<r- = i«-^a^)=c„e-^"^^-~S^^^ (2) 

where cjcj = or 1 is the fermion occupation number at 
site j. The expression for can be obtained by taking 
the Hermitian conjugate of . The string factor in the 
definition of a~ is necessary to ensure the correct anti- 
commutation relations between the fermionic operators, 
namely, {cm.cl} = Smn and {c„,c„} = 0. 

[In this paper, we will use both first quantized notation 
(wave functions) and second quantized notation (fermion 
creation and annihilation operators and occupation num- 
ber basis) as per convenience We specify the notation 
being used where necessary.] 

Following the Jordan- Wigner transformation, Eq. ^ 
takes the form 

N 

^ = - [2(4iCn+l +4+iC„) -I- 2hn cicn], (3) 

7l = -N 

where we have omitted a constant equal to . The 

fermionic operators can be represented in the momentum 
basis by the Fourier transform 

1 ^ 

71 — — N 

where the momentum k lies in the range [— tt, tt] and is 
quantized in units of 271 /{2N -|- 1); these operators sat- 
isfy the anti-commutation rules {cfc,c^,} = 6k,k' ■ In mo- 
mentum space, the first two terms of the Hamiltonian in 
Eq. ([3]) have the tight-binding form 

Hq = — ^ 4cos(fca) cJ,Cfc. (5) 

— 7r</e<7r 

As for the last term in Eq. ([3]), we consider the sit- 
uation where /i„ = hcos{Qn -\- 0). We note that this 
Hamiltonian appears in the Azbel-Hofstadter probleni^i 
of an electron hopping between the sites of a square lat- 
tice in the presence of a magnetic field applied in the 
perpendicular direction; the magnetic flux per square is 
proportional to Q. This problem has been studied ex- 
tensively, and we will refer to only some of the papers 
hero^^i^^. 

For the case Q = 7r/q, where q is an integer, the period- 
icity of the magnetic field is 2q. Using the decomposition 

we see that this term couples two fermionic modes with 
momenta fci and k2 if /ci = fc2 i 7r/g. This periodic term 
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FIG. 1: Picture of intermediate states for the case q — 2. 
The zero energy states at momenta k = ±7r/2 are connected 
to each other through the states lying at fc = and tt, as 
shown by red (dashed) and green (solid) lines respectively; 
the dispersion relation Ek ~ — cos k is indicated by a thick 
solid line. The system is described by a set of decoupled 
Hamiltonians Hk labeled by k lying in the range [— tt, —tt/2]; 
each Hk is a four-dimensional matrix in which the momenta 
k, k + 7r/2, k + TT and k + 3tt/2, lying in the regions separated 
by the vertical dotted hues, are coupled to each other. 

fragments the Hamiltonian into a 2q dimensional matrix 
form composed of momentum regions —Tr + rir/q < k < 
-TT + (r + l)7r/g with < r < 2g - 1. The tight-binding 
term of Eq. is diagonal in this basis while the periodic 
potential connects neighboring momentum regions. The 
corresponding matrix elements of Hk are given by 

(fc + rn/q\Hk\k + sir/q) 

= — 4cos(fc -I- sn/q) Sr.s - he^'^ 5r,s+i ~ /ie~"^ <5r.s-i, 

(7) 

where < r, s < 2q — 1, and we have assumed 'periodic 
boundary conditions' for the matrix Hk, so that r = — 1 
and 2q mean r = 2q—l and respectively. We have thus 
reduced Eq. ([3]) to the decoupled form 

— 7r+7r/(j 

H = ^ i?.. (8) 

k— — 7r 

Note that the total fermion number J\f — X)^iL-7v = 
c|.Cfc commutes with each of the Hk and is there- 
fore conserved in time. 

We now consider the energy spectrum of Hk, focus- 
ing on the low energy excitations that give the dominant 
contribution to the quench. As will be discussed in Sec. 
Ill, the initial condition for the quenching dynamics is 
such that the system is at half-filling at all times; we 
will therefore be particularly interested in the states near 
zero energy which correspond to the momenta k = ikp, 
where the Fermi momentum kp ~ tt/2. If the amplitude 
h of the magnetic field is zero, the system is gapless and 
the states at A: = ±7r/2 are degenerate with each other. 
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A small, non-vanishing value of h (compared with the 
bandwidth) breaks this degeneracy in a manner that we 
derive using a perturbative expansion in h. To consider 
coupling between the k = —tt/2 and fc = 7r/2 regions, one 
can go through one series of intermediate states lying at 
k ~ —TT/2+Tr/q, — 7r/2+27r/g, • • • , 7r/2— Tr/q (with an am- 
plitude equal to — /le"^ at each step), and through another 
series of intermediate states lying at fc = —tt/2 — n/q, 
— 7r/2 — 27r/q, • • • , 7r/2 + Tr/q (with an amplitude equal 
to — /ie~*'^ at each step). A picture of these two series of 
intermediate states is shown in Fig. [1] for q = 2. Each of 
these series consists of g — 1 intermediate states. At the 
q-th order in perturbation theory, we therefore obtain an 
effective Hamiltonian H^ff which has a matrix element 
between the states at fc = ±7r/2 given by 

A^{n/2\Heff\-n/2) = {-^ /2\H,f /2)* 

~ m=\ (4cos(-7r/2 + W9)) 
(-/ie-">)? 

^ mi! (4cos(-V2-Wz))' 

where the denominators come from factors like E_ 



(9) 



E 



-7r/2±s7r/g 



t/2 - 

4cos(— 7r/2 + ±s7r/(7) corresponding to the 



energies in the unperturbed Hamiltonian in Eq. ([5]). 
Simplifying Eq. © gives 

A = (V2|i/e//|-V2) = (-V2|ife//k/2)' 

g*90 _|_ ( — l)'J-lg-«90 



Yllll (4sin(^s/g)) 



(10) 



We thus find that the magnitude of A is given by 



h'i 2|cos(#)| 



4'"' n^;! sH^s/q) 

/i"? 2|sin((7<?!))| 
4^ Yllll Mns/q) 



if q is odd, 

if q is even. (11) 



In a sense, the phase (j) governs the relative phase between 
the two paths of intermediate states which connect the 
low lying states. We assume here that (j> is such that 
cos(g0) 7^ if g is odd, and sin(g0) 7^ if q is even; 
if these conditions are violated, we would have to go to 
higher order perturbation theory to find a non-zero ma- 
trix element connecting the states at fc = ±7r/2. 

We can now consider moving slightly away from fc — 
±7r/2; then the unperturbed energies of the states fc = 
—Tr/2 + k' and 7T/2 + k' arc given by — 4fc' and 4fc' respec- 
tively. The effective Hamiltonian describing these two 
states is then given by the 2x2 matrix 



H, 



eff,k' 



-4fc' A* 
A 4fc' 



(12) 



where we assume that A continues to be given by the 
expression in Eq. PH)) because fc' is small. The eigen- 
values of (IT^ arc given by ±y^ (4fc')2 + | Ap; this is the 



dispersion of a massive relativistic particle whose veloc- 
ity is equal to the Fermi velocity vp = 4, and mass is 
proportional to |A| ~ /i'' times cos{q(j)) or sm{q(f>). 

Hence, h = corresponds to a QCP where the mass 
gap vanishes. Given that the energy vanishes as |fc'| if 
h ~ and as h"^ if fc'' = 0, the dynamical critical exponent 
and correlation length exponent are given by z = 1 and 
V = q, respectively. The correlation length exponent v 
thus depends on the periodicity of the magnetic field 



III. QUENCHING DYNAMICS 

Having established the form of the Hamiltonian for the 
periodic system and its effective low-energy description, 
we now consider a specific quenching protocol. Given the 
magnetic field /i„ = hcos{Trn/q + (j)), we vary the am- 
plitude of the field in time as /i = t/r, where we refer 
to T as the quenching time. Att=— cx) (ft- = — 00), we 
start with the ground state of the Hamiltonian denoted 
by ^'o(— 00). (Throughout this section, the symbol ^ 
denotes first quantized wave functions). ^o(~c») is the 
state in which ct^ = 1 (cJjC„ = 1) at all the sites where 
cos{nn/q+(j)) > 0, and = —1 (c}jC„ = 0) at all the sites 
where cos{nn/q + 0) < 0. [If cos(7rn/g + cf)) ~ for some 
value of n, the ground state is not unique in the limit 
h = —00, since the states with c}jC„ = and 1 arc then 
degenerate. We therefore assume that cos{Tm/q + (p) ^ 
for all values of n and that a finite number of values of (jj 
are avoided. This is equivalent to the conditions imposed 
on cos{q<f>) or sin{q(l)) after Eq. (|lip .] It is clear that in the 
range 1 < n < 2q, cos{Trn/q + (f) is positive for half the 
sites and negative for the other half. Hence, the ground 
state is half-filled in terms of fermions. We can write 
the ground state as the product of ground states of the 
2g-dimensional Hamiltonians Hk{h = —00), 



*o(-oo) = (g)*o,fc(-oo), 



(13) 



where 4'o,fc(— 00) is a first quantized wave function which 
is obtained as follows. In the limit h — > —00, the on-site 
term is much larger than the hopping term in Eq. ([3]); 
hence the Hamiltonian Hk becomes independent of fc in 
this limit. This implies that since the ground state is half- 
filled in real space, the states corresponding to the wave 
function 'I'o,fc(^ —00) are also half-filled for each value 
of fc. Thus 'I'o,fc(^ —00) denotes the wave function 
corresponding to the state in which the q negative energy 
states of Hk{h = —00) arc occupied by fermions and the 
q positive energy states of Hk{h — —00) are empty. 

The system evolves dynamically according to the time- 
dependent Schrodinger equation. The state at time t is 
once again given by a product 



i*fe(0, 



(14) 



where the first quantized wave function ^'fe(t) is obtained 
by using H)^ (t) for the time evolution. In the limit t — > 00 
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{h — >• oo), we reach the state ^'(cio) = 04'fe(oo). Since 
the Hamiltonian changes time at a finite rate 1/r, we ex- 
pect that \l/(oo) will differ from the ground state of the 
final Hamiltonian, 5'o(oo) = {^^'o,fc(oo) except in the 
adiabatic limit. Om- goal is to determine how far ^'(oo) 
is from the final ground state ^o(oo) as a function of r, 
quantifiable by the total excitation probability P. [Note 
that the final ground state 5'o(oo) is also half-filled; in- 
deed this is the reason for choosing ft,„ to have an even 
period, 2q^ so that both the initial ground state and the 
final ground state have the same filling. If the period of 
hn were odd, the initial and final ground states, corre- 
sponding to /i = — oo and oo, respectively, would have 
different occupation numbers, and it would not be possi- 
ble for 5'o(— oo) to evolve in time to ^'o(oo) even in the 
limit r oo, since the dynamics conserves the fcrmion 
number] . 

The formal procedure for evaluating the final state, as 
employed in our numerical calculations is as follows. We 
define a 2(7-dimcnsional matrix M whose elements are 
given by 



Mrs 



Sr. 



s+1 



Sr,s-1 



(15) 



for < s < 2(7 — 1, and we again assume 'periodic bound- 
ary conditions' for the matrix M. In the limit h — oo, 
the magnetic term dominates the system and the eigen- 
values of the matrix Hk in Eq. (O are the same as those 
of M, while in the limit oi h oo, the eigenvalues of 
Hk are the same as those of —M. Let the q eigenvectors 
of M corresponding to the negative eigenvalues of M be 
denoted by ipi, where i = 1, 2, • • • ,5; we take the ij^i to 
form an orthonormal set. For /i — —00, the ground 
state is one in which the q states corresponding to the 
wave functions ■(/;; are occupied and the remaining q states 
(corresponding to the positive eigenvalues of M) are un- 
occupied. Next, let us assume that under time evolution 
from t = —00 to -l-oo using Hk{t), the q wave functions 
Tpi evolve into ■(/'i,fc(oo). Namely, 



■i/'j,fe(oo) = T exp -i 



(16) 



where T denotes the time-ordering symbol which is re- 
quired because Hk(t) varies in time since h = t/r. At 
t = 00, the ground state of Hk is one in which the 
states corresponding to the wave functions Tpi are excited 
states and are therefore unoccupied, while the states cor- 
responding to the remaining cigenstates of M are occu- 
pied. The probability of being in an excited state at 
i = 00 is then given by 



(17) 



(With this definition, pk always hes in the range [0,g]). 
Finally we obtain the total excitation probability P by 
integrating over k, namely. 



P = 



— Vk- 



(18) 



This is a measure of how far 'I'(oo) is from the ground 
state ^'o(oo) at i = 00; if P = 0, *(oo) = *o(oo) up to 
a phase. (We have defined the normalization in Eq. (fT8|) 
in such a way that P = 1 if = q for all values of 
W). In terms of the spins cr^ at different sites, P tells us 
the number of spins per site which point in the wrong 
direction at t = 00, i.e., in the direction opposite to that 
given by the energetically favorable state for Eq. (jS]) in 
the limit h = 00. Thus the total excitation probability P 
is related to the density of defects in the final state, i.e., 
density of spins pointing in the wrong direction. 

Analytical treatment:- We first capture the broad 
features of the quenching dynamics analytically by focus- 
ing on small h before delving into a more detailed numer- 
ical analysis. Although the quench tunes from h = —00 
to = 00, the excitations arc mainly produced during 
the time when h is close to zero, i.e., when \h\ is much 
smaller than the band width. This is because h ^ 
corresponds to a QCP where there are states lying ar- 
bitrarily close to zero energy. No matter how large the 
quenching time r, there are states whose energy is less 
than l/r; these are the states for which the dynamics is 
not adiabatic, and hence contributing significantly to the 
excitation probability. 

We thus consider the quenching problem in the basis 
of the effective Hamiltonian given in Eq. (|12|) for the two 
states at A; = —tt/2 + k' and 7r/2 -f A;' obtained by per- 
turbation for small h. As pointed out earlier, |A| scales 
as h'^ times cos{q(j}) or sa\{q(j)) depending on whether q is 
odd or even. Now, if h is varied in time as t/r, the time- 
dependent Schrodinger equation for the first quantized 
wave functions uik and U2k of each pair of momentum 
states takes the form 



dt 



Ulk 
U2k 



-4fc' (t/r)V(#) 
(t/r)«r(#) 4fc' 



Ulk 
U2k 



(19) 

where /(qtp) is equal to cos((/)) or sin(0) times factors 
given in Eq. (|11|) which are independent of r and 0. 
We can assume that f{q4>) is real and positive by gaug- 
ing away any phase dependence by a relative phase ro- 
tation between uik and U2k- Multiplying Eq. (jl9p by 
(^V/(#))^/^'^+^^ and defining t' = t{f {q<P) / r'tf / 
we obtain 



Ulk 

■ dt' \ U2k 

^ -4k' 



f{q<l>) 
t'l 



1/(9+1) 



t'l 



Ak 



' ( /(90) ) 



1/(9+1 




As a specific case, let us assume that q is odd; then 
t"^ — > ±00 for t' — > ±00 respectively. The ground 
state of the 2x2 Hamiltonian in Eq. ((20)) is then 
given by V- = (l/"\/2)(l, 1)"^ for t' —00 and by 
-0+ = (1/V2)(1,-1)^ for t' ^- 00 (here T denotes the 
transpose of a row vector). It is now clear that if we 
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start in the state ip- at t' = — oo, the excitation proba- 
bihty Pk = \{uik{t'),U2k{t'))*'4'+\'^ at t' = oo must be a 
function of a single variable given by A:'(r'//(q(/)))^/('+^^ 
Furthermore, by general quantum mechanical arguments, 
the excitation probability pk must be zero for r — >■ c» 
(adiabatic limit) and 1 for r (limit of sudden change 
in the Hamiltonian). Hence, if t is large, pk gains a 
significant contribution only from values of k' for which 

|fc'l(rV/(#))'/<''+'' < 1, i.e., Ifc'l < (/(#)/r'')^/('+^^ 
The exact form of pk depends on various parameters and 
has characteristic oscillations which are analyzed in the 
Appendix as well as in the numerics below. 

It now follows that the total excitation probability is 
given by 



P 



-7r/2+7r/(2g) 



Pk 



dk' 



k' 



(/M)) 



/(#) 

t1 



1/(9+1) 



(21) 



The second equation in Eq. ((2T|) follows if r is large 
enough so that the range of k' , which is of the order of 
(/('Z</')/'''')^^^'^^^'i is much less than the width of the mo- 
mentum band, ir/q; this justifies changing the limits of 
the integral in ^ from [-n /2 ~ t: / {2q), -tt /2 + n / {2q)] 
to [—00,00]. The third equation in (f2T|) follows because 
Pk is very small unless k' lies in a range of the order of 
ifiq(t>)/T''y^'''^~^^^- We thus see that t should be much 
larger than h/ J for the power-law scaling in the last equa- 
tion in dH]) to hold. 

Results:- Equation (PT|) is our central result. Con- 
tained in it is the l/r'^/'^'i+i) scaling behavior of the total 
excitations probability arising from spatial periodicity. 
The dependence on qcj) shows that the original off-set in 
the potential ft.„ = hcos{'nn/ q + </>) controls the amount 
of mixing between the various states, as argued in the 
perturbative derivation of the effective Hamiltonian of 
Eq. ((T^ . While this off-set does not contribute to the 
scaling form, it determines the magnitude of the proba- 
bility amplitude that shifts between the different states 
due to the quench. 

Numerical treatment:- We now test our results for 
various values of r, q and (p and study the detailed behav- 
ior of the quenching dynamics through numerical simula- 
tions. The procedure has been described above. Briefly, 
we begin with a state ipi corresponding to one of the neg- 
ative eigenvalues of the matrix M in Eq. ([T5|) . We then 
use Hk{t) given in Eq. (O, with h ~ t/T, to evolve ipi 
as shown in Eq. (jTB]) to find the state -0i_fc(oo). After 
repeating this calculation for all i from 1 to g, we calcu- 
late the excitation probability pk following Eq. ([TT]) . We 
then integrate over k as in Eq. ([18]) to obtain the total 
excitation probability P. 

One remark on the scaling forms for odd versus even 
values of q is in order here. As made explicit in Eq. ([5|), 




FIG. 2: Logarithmic plot of P versus sin{2(f>) for g = 2 and 
r = 8. The linear fit to the points has a slope of 0.31. 



the Hamiltonian Hk in Eq. ([7]) connects together all val- 
ues of k which are separated by multiples of 7r/g, and k 
lies in the range [— tt, — tt + tt /q]. This is the range of k 
considered for our numerical calculations. For the low- 
lying states, the value of k which lies in the above range 
and which is also connected to — 7r/2 is given by fco = — tt 
if q is even and fco = —tt + tt/ {2q) if q is odd. We therefore 
work with the scaling variable {k—ko){T'^/ sin(q(^))^/'^'"''^) 
if q is even and {k — ko){T'^ / cos{q(j))y^^'''^^^ if q is odd. 

For g = 1, we have /i„ = (t/T)cos{TTn). This prob- 
lem can be solved analytically as discussed in Ref. 0. 
For hn cos(7rn), there is a direct matrix element between 
states at momenta k and k-\-TT. Hence, for every value of 
k lying in the range [— tt, 0], we have a two-level system 
for which an exact expression for the excitation proba- 
bility can be found using the Landau-Zener treatment^^. 
We find that pk = e^^'^'^^™ ^ and the total excitation 
probability is 



P, 



q=l 



dk 

Pk 



(22) 



ff T > 1. 

For 9 = 2, the arguments above indicate that P should 
scale as {sm{2(/))y/^ /t^^^ . For the special case = 7r/4, 
the dependence of P on r was studied in Ref. [l3; nu- 
merical studies confirmed the —2/3 power law. We now 
present the dependence of P on sin(2(^) in Fig. [2 a lin- 
ear fit to the logarithmic plot gives a slope of 0.31 which 
is close to a 1/3 power law. The fact that P exhibits 
a power-law scaling with respect to both r and sin(20) 
(for q = 2) confirms the validity of the perturbation the- 
ory developed in Eqs. ([ WTI]) for h close to zero and the 
fact that the excitation probability is dominated by the 
behavior for small h. 

For q = 3, our arguments show that P should scale 
as {cos{3(l>)y^'^/T^^^. Fig. |3] shows a logarithmic plot of 
P versus r for t/i = 0; a linear fit gives a slope of —0.77 
which is in good agreement with the expected value of 
—3/4. Figure H] shows a plot of pk versus {k + 57r/6)r'^/^ 
for three values of r, for </> = 0; the three curves are seen 
to coincide, indicating that pk is indeed a function of the 
scaling variable (fc -t- 5tt/6)t^^^ as indicated in Eq. (|20p . 
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FIG. 3: Logarithmic plot of P versus r for q — 3 and = 0. 
The linear fit to the points has a slope of —0.77. 
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FIG. 4: Plot of pfe versus (fc + 5n/6)T^''^ for r = 2 (red), 4 
(blue) and 8 (black), for q — 3 and (f> = 0. 



and only a range of values of fc -|~ 57r/6 of the order of 
1/^3/4 -g gggj^ contribute significantly to the total ex- 
citation probability. 

Figure |4] shows several oscillations in as a function 
of k. To obtain a clearer understanding of these oscil- 
lations, we study in the Appendix a two-level system in 
which the one of the terms in the Hamiltonians varies 
in time as \t\'^sgn{t); this is motivated by the form of 
the Hamiltonian in Eq. ((20)) . We find that the number 
of oscillations increases with q. Our analysis also pro- 
vides a derivation of the period of oscillation for large q, 
indicating that the oscillations in k have a period that 
also scales as l/r'?/('?+i)- ^j^jg follows from the statement 
proved in the Appendix that, for large values of q, the 
oscillations have a period 7r/2 in a parameter b which is 
equal to the quantity 4fc'(rV/(#))^/^'+^^ in Eq. 
To see the manner in which pk depends on q, wc have 
plotted Pk versus k + 2l7r/22 for q = 11 and </> = in 
Fig. [5] Note that for g = 11, Pfc is expected to be a func- 



1.4p 
1.2- 



1 - 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 

k + 21ji/22 

FIG. 5; Plot of pfe versus k + 2l7r/22 for q = 11, (/> = 0, and 
T = 16. 

tion of (fc-K2l7r/22)T"/i2; the fi gure indicates that only a 
range of values of fc-|-2l7r/22 of the order of I/t^-^^^^ con- 
tributes substantially to the excitation probability. The 
number of oscillations in pk shows a clear increase com- 
pared to the case of g = 3 in Fig. S) 

Our numerical simulations thus confirm the scaling be- 
havior of P as a function of the quenching rate and off-set 
0, showing that the scaling holds for a significant range 
of parameter space. They also illustrate the exact man- 
ner in which p^ varies as a function of momentum k and 
periodicity q. 



IV. MORE COMPLEX PERIODIC 
STRUCTURES 

In the previous sections, we have considered a magnetic 
field given by /i„ = hcos(TTn/q+(f>) assuming even period- 
icity 2q, and we have varied h in time as t/r. Our analysis 
of the energy spectrum close to the points fc = ±7r/2 and 
the scaling of the excitation probability were based on 
perturbation theory. We can now ask if other periodic 
forms for hn lead to similar results. 

We again assume that /i„ = hcos{Qn + cj)), but we set 
Q ~ Trp/g, where Q < 27r, p and q are relatively prime 
integers, and p is odd. These conditions ensure that the 
period of hn is an even integer given by 2g, so that the 
ground states for h very large and either positive or neg- 
ative are both half-filled. Once again we can do pertur- 
bation theory to obtain the matrix element between the 
states at fc = ±7r/2; we find that at the lowest order in 
h which is given by h"^, the magnitude of the matrix ele- 
ment is given by Eq. (jlip regardless of the value of p. This 
follows from the number theoretic fact that the sets of in- 
tegers (1, 2, • • • , q — 1) and {p, 2p, • • • ,p{q — 1)) modulo q 
are identical if p and q are relatively primo^l; thus the set 
of positive numbers | sm{T:sp/q)\, for s = 1, 2, • • • , q~l, is 
the same for all values of p which are relatively prime to 
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FIG. 6: Plot of pk versus k + n for 
r = 8. 
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FIG. 7: Plot of Pk versus k + n for 

T = 8. 



q. Hence the effective two-level problem given in Eq. ([19]) 
depends only on q and not on p. 

As a simple example, we compare the cases given by 
Q = it/4: and Q = 3tt/A. For both of these, /i„ has 
period 8, and for small h, perturbation theory leads to 
the same expression for the magnitude of the matrix ele- 
ment given in Eq. pT|) between the states at fc = ±7r/2, 
namely, |A| = (l/16)/i^| sin(4(/))|; we therefore obtain the 
same effective two- level problem given in Eq. ([T5)) . How- 
ever, the complete 8-dimensional Hamiltonians Hk are 
not identical for Q — 7r/4 and 37r/4 for arbitrary values 
of h and fc, even if we allow for unitary transformations. 
It is therefore interesting to compare the results for the 
excitation probabilities pk in these two cases. This is 
shown in Figs. [H] and [TJ we see that the forms of pk are 
qualitatively very similar although they differ quantita- 
tively in the two cases. 

Next, we consider the situation in which the magnetic 



FIG. 8: Logarithmic plot of P versus r for hn = 
(t/r)[0.2cos(7rn) -|- cos(7rn/2 -|- 7r/4)] showing a cross-over be- 
tween two power laws. The two linear fits near the beginning 
and end of the plot have slopes of —0.65 (blue solid) and —0.48 
(red dashed) respectively. 



field hn is periodic in n, but has several Fourier compo- 
nents. (We continue to assume that the period of hn is 
even, with /i„ being positive on half the sites and neg- 
ative on the other half). For instance, let us consider a 
field of the form /i„ — h[ai cos(7rn)-|-a2COs(7rn/2-|-7r/4)], 
where ai, 02 and (f) are constants; we then vary h as t/r. 
(In order to to focus on the scaling with respect to the 
quenching rate, we have chosen the phase of the period 
4 term to be 7r/4 such that sm{2(j>) = 1). If the a^'s are 
small, the perturbation theory described in Sec. II shows 
that there is a coupling between the low-energy states 
at fc = ±7r/2 at first order in ai and at second order 
in 02; these would lead to a total excitation probability 
scaling as a}^^/r^/^ and (a-2)'^l^ jr^l^ respectively. As- 
suming that a\ and a-i are both non-zero, we expect that 
for very large values of r, the l/r^^^ scaling would dom- 
inate over the l/r^^'^ scaling. But if t is not too large, 
then the scaling could dominate if the values of 

T, fli and 02 are such that (a2)^/^/r^/^ > a-^"^ jr^l'^. In 
particular, if (02)^^'^ 3> aj^^, wc may expect to sec a 
1 / T^l'^ scaling over a range of values of r before crossing 
over to a scaling for very large values of r. Fig- 

ure [8] shows a logarithmic plot of P versus r for the case 
hn = (t/T)[0.2cos(7rn)-|-cos(7r7i/2 + 7r/4)] (corresponding 
to a\ = 0.2 and = 1), along with two linear fits near 
the beginning and end of the plot. The linear fits have 
slopes of —0.65 and —0.48 which are in good agreement 
with the values —2/3 and —1/2 respectively. 

Finally, we may ask what happens ii Q = Trp/q, where 
q is very large and the variation of the potential is very 
slow. (For instance, this would be useful if we were in- 
terested in studying the case of a quasiperiodic pattern 
of the field hn, and we approached such a pattern by 
taking rational approximations of the form p/q with q 
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becoming progressively large). Unfortunately, it seems 
to be rather difficult to study the case of very large q ei- 
ther analytically or numerically. Perturbation theory up 
to g-th order in h would not be reliable, since the ener- 
gies appearing in the denominator of Eq. (jlO[) . sm{T: s/q), 
come arbitrarily close to zero if s is held fixed and g — > oo. 
Numerical calculations based on time evolution with the 
2g-dimensional Hamiltonian Hk in Eq. ([7]), with h = t/r, 
also become difficult if q is very large since there would be 
a large number of levels which come close to each other at 
different times; the time evolution would therefore need 
to be done very accurately to go from t = — oo to t = oo. 
In this limit, further studies and alternate methods are 
called for. 



V. QUENCHING IN A 
TOMONAGA-LUTTINGER LIQUID 

We have studied quenching in a periodic system of 
free fcrmions; we now turn to the manner in which 
this physics becomes modified in the presence of interac- 
tions between fermions. Our arguments borrow from the 
quenching treatment of Ref. |lj for interacting bosons 
in a periodic potential favoring one particle per poten- 
tial minimum, adapting it to the situation of interacting 
fermions as well as generalizing to multiple boson occu- 
pancy. 

Tracing back to our initial starting point, the spin 
Hamiltonian of Eq. ([Ij, the addition of a nearest- 
neighbor zz coupling of the form JzC^a^^i results in an 
interaction term in the fermionic system. Employing the 
Jordan- Wigner transformation described in Sec. II leads 
to the addition of a four-fermion interaction of the form 



in the transformed fermionic Hamiltonian 



of Eq. ©. 

We can treat the effect of these interactions on the 
quenching process described in the previous sections as 
follows. The quenching involves tuning the strength of 
the periodic potential h from — oo to +00 as a function of 
time, and once again, the most important contribution 
to quenching dynamics comes from the gapless region at 
vanishing h. We can thus consider the low-energy physics 
of the fermionic system for small h, namely the effective 
Hamiltonian Eq. (|12p which consists of right-moving and 
left-moving linear modes at Fermi points k = ±tt/2, t/ijj 
and V'L J respectively, and a coupling between them of or- 
der /i' derived from q-th order perturbation. In addition, 
we now have interactions and these can be treated within 
the standard context of bosonization, where the right and 
left moving modes V'fl/L ^ exp[i^ (=F\/X$ + 6/VK)] 
are expressed in terms of the bosonic fields $ and 6^. 
In terms of the scalar field $, the deviation of the den- 
sity of fermions p from the average density po is given 
by P — Po = —{\/K/'K)d^/dx. The Luttinger param- 
eter K characterizes the strength of the interactions 
{K < (>)1 for repulsive (attractive) interactions, K = 1 
in the absence of interactions.) For a mapping from 



the spin model, the procedure can be considered in the 
range — 1 < Jz ^ 1; the parameter K is then given by 
so that K goes from 00 to 1/2 as 



K — 2 

~ 2cos-i(-,/^) ' 

goes from —1 to 1. 

The linear modes at /c = ±7r/2 and the short-range 
fermion interactions together provide the standard action 
describing a Tomonaga-Luttinger liquid (TLL): 



S = 



dtdx 



'dt 



9$ 

dx 



(23) 



(Here, the velocity of the bosonic excitations has been 
set to unity.) The contribution of the quench term to 
the action has the sine- Gordon form 



Srr 



j j dtdx h'^ cos(2V7r¥<I>), 



(24) 



Compared to the standard cases of interacting fermions 
that result in forms similar to Eq. and (P^ . we em- 
phasize that these terms have been derived as the low- 
energy description of a system having a higher periodicity 
2(7, resulting in the key dependence of the mass term Sm 
on q. 

Using Eqs. ((23)) and (|24)) as the starting point, we 
invoke renormalization group (RG) arguments^^ to de- 
rive the scaling of defect formation due to the quench. 
By computing the correlation function of the opera- 
tor cos(2V7r i4r$) at two different space-time points, one 
can show that it has mass dimension K; let us denote 
the coefficient of this operator in the action by p, in 
general. The parameter p effectively becomes depen- 
dent on the length scale L, satisfying the RG equation 
dp/dlnL = {2—K)p. Given the initial value of p{a) = /i' 
at some microscopic length scale a (this could be either 
the lattice spacing or the average distance between the 
particles in a continuum model), the solution of the RG 
equation is p{L) = h'^{L / af'^^ . Assuming that h'' <^1, 
we see that p{L) is of order 1 at a length scale given 
by L = ^, where {i/a) - l/Zi'/^^-if)^ xhe condition 
that p{^) ^ 1 implies that we have reached a strong cou- 
pling (disordered) regime at the length scale ^, i.e., ^ is 
the correlation length of the theory. Since h denotes the 
deviation from the quantum critical point, the relation 
^/a ~ l/Zi'/'^^"^) implies that the correlation length ex- 
ponent is given by v = ql{2 — K). As argued in the 
previous sections, if h is quenched through the QCP at 
h — at a, rate given by 1 /r, the KZ scaling form then im- 
plies that the density of defects will scale as 1/t'^'^/'-^'^'^^\ 
Hence, substituting the value of jy in this scaling form, 
the density of defects scales as 



PD ^ 1/t 



q/{q+2-K) 



(25) 



This scaling behavior can be derived more rig orously by 
trivially modifying the treatment in Ref. [ij to include 
the q dependence. For the non-interacting case, K = 1, 
this form reproduces the earlier expression 
In the presence of interactions, we see that the power 
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varies continuously with K. We note that the power law 
ll^q/(q+2-K) Qj^jy ^^Yid for 1/2 < if < 2. The value 
K ~ 2 corresponds to a Kosterlitz-Thouless transition. 
For if > 2, the cosine term in the action is irrelevant. 
It has been argued in Ref. [l2| that the probabihty of 
excitations then receives contributions from all modes, 
not just the low-energy modes near k = ±7r/2; hence the 
scaling law has the form l/r regardless of the value of q. 

These arg uments can be used to generalize the results 
of Ref. [ij to a multi-boson situation. In addition to 
fcrmionic systems, the TLL form can be equally well ap- 
plied to one-dimensional systems of repulsively interact- 
ing bosons; the Luttingcr parameter K goes from oo to 1 
as the strength of the repulsive interaction is varied from 
zero to oo. In particular, when the interaction strength is 
infinitely large, the bosonic system is equivalent to a sys- 
tem of non-interacting spinless fermions^^ with K = 1. 
In the continuum of interacting bosons described by the 
action in Eq. (|23p , one can introduce a spatially periodic 
potential of the form V{x,t) = h{t) cos{2tt px / q) , where 
q is an integer. For very large values of h{t), the ground 
state would have exactly q particles residing in each of 
the minima of this potential. (Therein lies the general- 
ization of Ref. [ill which assumes q = 1.) To analyze a 
quench of the form h ~ t/r, we note that there is a one- 
to-one correspondence between the interacting system of 
bosons and that of the fermions studied above. Focusing 
on the lower energy physics, once again Eq. (j23p describes 
a gaplcss TLL with a Fermi momentum given hy kp = np 
while Eq. (|24p describes the contribution coming from a 
small, non- vanishing h. The subsequent analysis above 
also goes through for the bosonic case except for a redef- 
inition of the interaction parameter. We thus conclude 
that the density of defects (which corresponds to some of 
the potential wells having either more than or less than q 
particles) scales as i/r'^/ii+^-K) ^ ^ 2, 

and as l/r if K > 2. 



VI. CONCLUSION 

To summarize, we have studied the dynamics of a spin- 
1/2 chain in which the amplitude of a spatially periodic 
magnetic field h cos (Qn -1-0) is slowly varied in time so as 
to take the system across a quantum critical point. We 
have equivalently analyzed a system of spinless fermions 
with a spatially periodic and time-dependent chemical 
potential. We find that a quenching rate of l/r takes the 
system to a final state having a density of defects that 
scales as 1/t'^/'^i+^^ for Q = 7r/g; the power depends on 
the period of the magnetic field. We have shown this by 
using perturbation theory to derive an effective Hamilto- 
nian for pairs of states near zero energy; this Hamiltonian 
has a parameter varying as a power of the time, where 
this power also depends on the period of the field. We 
have confirmed our results numerically in a number of 
cases. If the magnetic field has several Fourier compo- 
nents, the power law corresponding to the smallest pe- 



riod is found to dominate for very small values of 1/t 
although one may find a cross-over at intermediate val- 
ues of T depending on the amplitudes of the different 
Fourier components. If the spin Hamiltonian has addi- 
tional terms which can be mapped to interactions be- 
tween the fermions, the power varies continuously with 
the interaction strength. To obtain a better understand- 
ing of some features of the excitation probability in the 
spin chain problem, we have numerically studied a two- 
level system in which a term in the Hamiltonian varies as 
a power law in time. We find that the excitation proba- 
bility in this problem has a complicated dependence on 
the power resembling the results obtained for the many- 
body spin or fermion systems. 

Our predictions, of interest from the perspective of 
critical phenomena and out-of-equilibrium dynamics, can 
be investigated in various physical systems. The most 
immediate experimental realizations would perhaps be in 
the area of cold atoms or molecules trapped in effectively 
one-dimensional optical lattices^, where tuning parame- 
ters is easy and quenching has become an extensive topic 
of study. Another area would be one that has been stud- 
ied in great depth and well characterized, namely, various 
physical systems described by spin chain physics. In this 
case, the only additional ingredient necessary would be 
to subject the chains to a magnetic field having a peri- 
odic modulation in space and whose strength could be 
dynamically tuned. As yet another instance, recent pro- 
posals have employed effectively spinless p-wave paired 
fermionic superconducting wires in schemes for topologi- 
cal quantum computation requiring tuning between topo- 
logical and non-topological regions in these wires. As 
described in our previous workii, periodic potentials in 
such systems offer a way of studying topological aspects. 
A systematic study of the quench work presented here 
generalized to include anomalous terms due to supercon- 
ductivity (already initiated in Ref. [13) would be useful 
for the proposed schemes and as a study in of itself on 
quenches in topological systems. 

Our results have revealed how mode coupling can re- 
tain the slow out-of-cquilibrium dynamics expected of 
quenching through a QCP while giving it a richer, more 
complex structure. Here, the mode coupling was via frag- 
mentation of the Brillouin zone by a periodic potential. 
Future avenues would involve more intricate mode cou- 
pling, for instance, in the presence of quasiperiodic poten- 
tials, or in disordered systems, which completely break 
translation invariance and can undergo dclocalization- 
localization transitions, or even in systems of weakly cou- 
pled isolated quantum states. Finally, our studies have all 
been in the zero temperature limit. Since the quantum 
critical regime extends up to some finite range of tem- 
peratures, we expect the results presented here to persist 
at temperatures which are much smaller than the energy 
gap at the initial timeii; a rigorous finite-temperature 
study is in order. 
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Appendix A: Non-linear quenching in the 
Landau-Zener system 

In this Appendix, we will study a generalization of 
the Landau-Zener proble m^^'^^ in which the Hamiltonian 
varies non-linearly in time. We consider a two-level sys- 
tem evolving with a time-dependent Hamiltonian as 



dt 



Ul 
U2 



\t\isgn{t) 
b 



Ul 
U2 



(Al) 



where 6 is a constant, and q is a positive integer. It 
is sufficient to study this problem for the case that b is 
real and positive, since the Hamiltonian for any complex 
value of b is related by a unitary transformation to a 
Hamiltonian in which b is replaced by |&|. The reason for 
studying the Hamiltonian in Eq. (|A1[) is that it is related 
by a unitary transformation to the one in Eq. ()20|) . with 
b playing the role of 2k'{Ti / f{q(t>)y/'^''+^\ Note that we 
have taken the time dependence in (|A1[) to be of the form 
\t\''sgn{t); this has been done so that the ground states 
for t — >■ —00 and t ^ oo are given by (ui, U2) = (1, 0) and 
(0, 1) respectively, just as they are for the case of linear 
time variation (q = 1). 

Note that we have taken the time dependence in (jAl[) 
to be of the form \t\''sgn{t), rather than the two ex- 
pressions agree for q odd, but not for q even. We have 
chosen the form \t\'^sgn{t) so that for all values of q, the 
ground states for t — > —00 and < — > 00 are given by 
(wi,W2) = (1,0) and (0,1) respectively, just as they are 
for the case of linear time variation (g = 1). 

At t = —CG, we begin with the ground state (mi, W2) — 
(1,0). We then numerically evolve the wave function as 
in Eq. (jAip up to i = 00; at that point, we find the ex- 
citation probability given by p{b) = |wi(oo)p. For 6 = 0, 
we expect p = 1 since the Hamiltonian has no matrix ele- 
ment connecting the initial and final ground states; hence 
the system remains in the initial ground state (the sudden 
limit). For & — > cxj, we expect p ~ since the instanta- 
neous ground and excited states remain well separated 
at all times; hence the system follows the instantaneous 
ground state (the adiabatic limit). For q = 1, one has 
the exact expression p{b) = exp(— 7r6^)^. For q > 2, the 
exact expression for p{b) is not known. 

Figure [S] shows how p varies with b for three different 
values of q. We see that p{b) is a product of an oscillatory 
function and a decaying function of b (this seems to be 





2 4 6 



FIG. 9: Plot of p versus 6 for g = 5 (top), q — 10 (middle), 
and q — 20 (bottom). 



true for all values oi q > 2). The period of oscillations 
appears to approach 7r/2 for large values of q and small 
values of b. We can understand this as follows. If q ^ 1, 
the function is much larger than 1 for |i| > 1 and 
much smaller than 1 for \t\ < 1. Let us therefore make 
the approximation — 00 for \t\ > 1 and = for \t\ < 1. 
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FIG. 10: Plot of p versus t iov q = 20 and 6=1. 
The Hamiltonian in Eq. ()Aip can then be written as 

Since we begin with the wave function {ui,U2) = (IjO) 
at i = — CX3, and the energy levels are given by ±oo from 
t ~ — oo up to t = 1, the wave function remains equal 
to (1,0) (times a phase factor due to time evolution with 
infinite energy) up to t = — 1. Then the Hamiltonian 
is independent of time from t = — ltot = l; this is 
known as the waiting problem^. It is easy to solve for the 
time evolution during this period; we find that the wave 
function at t = 1 is given by (cos(26), i sin(26)) times 
a phase factor. From t ~ —1 to t = oo, the energy 
levels are again given by ±oo; hence the wave function 
at t = oo is given by ui = cos(26) times a phase factor 
and U2 = sin(26)) times another phase factor. We thus 
obtain p{b) = cos^(26). This explains the oscillations in 
p versus b with period 7r/2 and amphtude equal to 1/2 
for large values of q and small values of b. However, 
this simple explanation breaks down for large values of 



b where the period of oscillations starts deviating from 
7r/2 and the amplitude of oscillations goes to zero as we 
can see in Fig. [9l 

Figure [To] confirms the scenario presented above for the 
variation of the excitation probability p = versus 
the time t. Namely, p does not change much before t = 
— 1 or after t = 1, but it oscillates between t = —1 and 1. 
Figure [11] shows a two-parameter fit to a plot of p versus b 
for g = 10; a function of the form p{b) = e"'^'^^ cos^(2c2&) 
with ci = 0.091 and C2 = 1.11 gives a good fit for b < 3.5, 
but the fit becomes progressively worse for larger values 
of b. For q = 20, we find that the same function of p 
versus b gives a good fit with ci = 0.020 and C2 = 1-08 




b 



FIG. 11: Two-parameter fit to plot of p versus b for q — 10. 
The fitting function p{b) = e'"^''^ cos^(2c2b) gives ci = 0.091 
and C2 = 1.11. 



for b < 5.5. Note that our analysis above predicts that 
C2 should be equal to 1. 

The study of the two-level system in this Appendix 
sheds some light on the results we obtained in Sec. HI for 
the spin system. In particular, the observation in Fig. [5] 
that as q increases, p{b) shows increased oscillation as 
a function of b (while the amplitude of the oscillations 
goes to zero for large b) explains similar features seen in 
Figs. H] and [3 for the spin system. 
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